Strong Electron-Phonon Coupling in Superconducting MgB2: A Specific Heat Study 
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We report on measurements of the specific heat of the recently discovered superconductor MgB2 in 
the temperature range between 3 and 220 K. Based on a modified Debye-Einstein model, we have 
achieved a rather accurate account of the lattice contribution to the specific heat, which allows us to 
separate the electronic contribution from the total measured specific heat. From our result for the 
electronic specific heat, we estimate the electron-phonon coupling constant A to be of the order of 
2, significantly enhanced compared to common weak-coupling values < 0.4. Our data also indicate 
that the electronic specific heat in the superconducting state of MgB2 can be accounted for by a 
conventional, s-wave type BCS-model. 

PACS numbers: 74.25.Bt, 74.60.-w, 74.25.Kc 



The recent discovery of superconductivity in MgB2 be- 
low Tc « 39 K j^] has caused a remarkable excitement in 
the solid-state physics community. Critical temperatures 
of this magnitude inevitably raise the question whether 
mechanisms other than the common electron-phonon in- 
teraction are responsible for the transition. In their very 
recent work, Bud'ko et al. ||] have investigated the Boron 
isotope effect in superconducting MgB2 and found that 
replacing ^^B by ^°B increases the critical temperature 
by about 1 K. This was taken as strong evidence that su- 
perconductivity in MgB2 is of conventional nature, i.e., 
the electron pairing interaction is phonon-mediated. Ko- 
tagawa et al. have interpreted their ^^B NMR mea- 
surements as indicating strong coupling s-wave super- 
conductivity 1^. Preliminary ^^B NMR measurements 
at low temperatures yield spectra which are consistent 
with the expectations for a type-II superconductor in the 
mixed state Q. Evidence for sizeable electron-phonon 
coupling is also provided by recent band-structure calcu- 
lations 1^-01 . Various tunneling experiments have pro- 
vided some evidence for conventional BCS-like supercon- 
ductivity, but the values of the superconducting energy 
gap extracted from these measurements vary from 2 to 7 
meVH^Tg. 

In this letter, we report on measurements of the spe- 
cific heat Cp of MgB2 and present Cp{T) data in a wide 
temperature range. Using a consistent model for the con- 
tribution of the lattice vibrations to Cp{T) we calculate 
the electronic specific heat of this material. We show that 
the electronic specific heat below Tc is in good agreement 
with a conventional BCS-type interpretation. 

The sample has been prepared from commercially 
available MgB2-powder (Alfa Aesar) by sintering a 
pressed pellet at 500°C for about 72 h. Electron mi- 
croprobe investigations of the sample have shown that 
impurities of heavy elements (Cu, Ni, W) with concen- 
trations of the order of 10~^ are present in the sample. 

The specific heat Cp{T) of the sintered MgB2 sample 
has been measured using two different experimental tech- 
niques in overlapping temperature ranges. A standard 
relaxation technique was employed in the temperature 



range between 3 and 45 K. For temperatures between 
20 and 220 K an adiabatic continuous heating calorime- 
ter was used. Special care was taken to minimize the 
radiation losses at elevated temperatures The tem- 



peratures in the range between 3 and 45 K were reached 
using a pumped "^He cryostat, and for those between 20 
and 220 K, a conventional gas-fiow '^He cryostat was used. 

In the lower inset of Fig. |^ we show the magnetisa- 
tion of our sample divided by the constant applied field, 
M/H, as a function of temperature for T < 50 K. The 
superconducting transition temperature is 37.5 K and is 
indicated by the vertical arrow in the figure. The reduc- 
tion of Tc of our sample, compared to Tc « 39 K observed 
by Bud'ko and coworkers |^ , is most likely due to the fact 
that our sample is not as clean as theirs and the super- 
conducting transition temperature is slightly reduced by 
the impurities. We note, however, that the anomaly in 
Cp(T) at Tc is at least as sharp as the one reported in 
Ref. [|. 

In Fig. |l| we show the as measured specific heat Cp (T) 
vs. T in the whole temperature range covered in this 
study. The upper inset of this figure shows the same data 
in a limited temperature range around Tc- The anomaly 
in the specific heat due to the superconducting transition 
of MgB2 at about 37.5 K is clearly resolved. We note that 
the absolute values of the data presented here arc in good 
agreement with previously reported results of the specific 
heat of MgB2 in a narrow range of temperatures around 

Tc i- 

With respect to superconductivity in MgB2, the most 
interesting part of Cp{T) is the electronic specific heat 
Cei{T). In order to reliably separate this contribution 
from the total, measured specific heat, the contribution 
due to lattice excitations, Cph (T) , has to be known quite 
accurately. As may be seen in Fig. ^ the lattice provides 
the dominant contribution to the total, measured specific 
heat above 20 K. Well below the Debye-temperature On, 
the specific heat in the normal state of a common metal 
is usually approximated by CpiT) = ^T + f3T'^, where 
the first term represents the electronic and the second 
term the lattice specific heat. We show below that this 
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approximation for evaluating both, the electronic and the 
lattice specific heat, is not applicable for MgB2 at T > Tc- 
Even at temperatures only slightly above Tc, the lattice 
specific heat may not simply be described by assuming a 
linear dispersion of the acoustic modes. 
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FIG. 1. Specific heat Cp{T) of sintered MgB2 as a function 
of temperature between 3 and 220 K. The upper inset shows 
the same quantity in a limited temperature range. The lower 
inset displays the magnetisation M divided by the applied 
field // as a function of temperature for our sample. The 
vertical arrow marks the onset of superconductivity in this 
material at Tc ~ 37.5 K. 



In Fig. 1^, we show the total measured specific heat of 
MgB2 divided by as a function of T. We note that, 
just above Tc, Cp(T)/T^ increases with increasing tem- 
perature and passes through a pronounced local maxi- 
mum at about 60 K. Such a feature cannot be described 
by the approximation Cp(T) = 7T + PT^ , mentioned 
above. 

Due to the present lack of thermal expansion data, 
we cannot calculate the specific heat Cv(T) at constant 
volume. The difference between Cp{T) and Cy(T) is 
expected to be small in the entire covered temperature 
range and at this point we neglect it. 

We assume that the total specific heat at T > Tc may 
be described by an electronic contribution Cei(T) which 
is given by 7T at T > Tc in the whole covered tempera- 
ture range, and a lattice contribution Cph(T). 

The lattice specific heat is generally given by 
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where g{iv) denotes the phonon density of states (PDOS), 
h the Planck and kg the Boltzmann constant, respec- 
tively. In the Debye approximation, the lattice is treated 
as an isotropic continuum with a linear dispersion, which 
leads to a PDOS proportional to for uj < ujd, ^^d 
denoting the cut-off frequency above which the PDOS is 



zero, and accordingly, to a low-temperature lattice spe- 
cific heat of the form Cph{T)/T^ =constant. As men- 
tioned above, it is not possible to describe the observed 
maximum in Cp{T)/T^ by this simple model. The Debye 
model may easily be extended to include deviations from 
the linear dispersion. This leads to a PDOS of the form 
17(0;) — fiuj^ -|- vuj'^ and, in the low temperature limit, to 
Cph/T^ = /3 + 5T^ . The cut-off frequency lou is chosen 
such that the total PDOS per mole is limited to 2>pNA^ 
where p is the number of atoms per unit cell and Na 
Avogadro's constant. The lattice specific heat in this ex- 
tended Debye-scheme is then given by 
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Kortus and coworkers |^ have recently computed the 
energies of the zone-center optical modes in MgB2 using 
a frozen phonon scheme. According to their work, the 
lowest optical mode is doubly degenerate and located 
around an energy of hajopt/ks ~ 460 K. A calculation 
of the phonon spectrum yields a peak in the PDOS at 
Tiujopt/kB ~ 380 K j6j. Modes of this energy contribute 
to the lattice specific heat already at temperatures far 
below Tc and thus cannot be neglected in our analysis. 
In order to take them into account, we include their con- 
tribution into our approximation of the lattice specific 
heat by treating them as Einstein modes. The specific 
heat of one Einstein-mode [Na states at an energy Tioje) 
is given by 



Copt{T) = Na 



exp(e£/T) 



ksT^ (exp(e£;/T) - 1)' 
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where = denotes the Einstein-temperature, and 
Tilje the energy of the optical mode. 

MgB2 has p — "i atoms per unit cell, and therefore 
9 phonon-modes. We treat the energetically lowest op- 
tical mode, which is doubly degenerate, according to 
Eq. (^) and the remaining 7 modes by the formula given 
in Eq. (||) . The solid line in Fig. ^ is a fit based on 

Cp{T)/T^ ^-ilT'' + Cd{T)/T^ + CoptiT)/T^ (4) 

to the data at T > Tc. By inspecting Fig. ^, it may 
be seen that with this model the measured total specific 
heat Cp{T) may be well approximated at T > Tc, thus 
providing a rather accurate description of the lattice and 
the electronic specific heat of MgB2 in its normal state. 

The fitting procedure provides a 
value of 5.5 mJ/mole • for 7. This value is almost 
twice as large as the 7 value given in Ref. However, 
considering the high accuracy of our fit we are confident 
that our value for 7 is quite reliable. If we interpret the 
7 parameter as the Sommerfeld constant, a comparison 



2 



50 100 150 200 

T(K) 

FIG. 2. Specific lieat Cp{T) of MgB2 plotted as Cp{T)/T^. 
The solid line represents a fit of Eq. (M) to the data at T > Tc, 
and the broken line represents the resulting lattice contribu- 
tion (see text). 
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FIG. 3. Total specific heat Cp{T) minus the lattice con- 
tibution Cph(T) divided by T vs. temperature. The solid 
line represents a fit to our data using Eq. The broken 

line represents the resonance contribution due to the impuri- 
ties and the dash-dotted line the calculated electronic specific 
heat (see text). 



with recently calculated densities of states at the Fermi- 
energy Ef, D{Ef) = 0.72 states/unit-cell-eV g and 
0.74 states/unit-cell-eV ||l2[, respectively, leads directly 
to a value ^ = 3.14 for the average mass-enhancement 
of the conduction electrons in MgB2. Neglecting other 
many-body effects, the corresponding electron-phonon 

coupling constant A = ^ 1 = 2.14 is significantly 

enhanced above the usual BCS- weak-coupling values of 
Xwc < 0-4, thus providing evidence for a considerable 
electron-phonon coupling in MgB2. 

The lattice contribution to the specific heat extracted 
from the fit is shown in Fig. ^ as a broken line. Since 
our model reproduces the total specific heat at elevated 
temperatures rather well, we may safely assume that our 
calculation provides the lattice contribution not only for 
T > Tc, but across the whole covered temperature range. 
The other two main parameters which emerge from the 
fit are the Debye temperature 9d = 746 K in good 
agreement with the value given in Ref. §, and the Ein- 
stein temperature 9e — 325 K. The corresponding energy 
for this dispersionless mode is in reasonable agreement 
with the position of the peak in the PDOS calculated in 
Ref. and observed at huj/kB — 365 K in the inelastic 
neutron scattering data of Ref. . 

With the proviso that the lattice specific heat Cph{T) 
is now established at all covered temperatures, we calcu- 
lated the electronic specific heat of MgB2 by subtracting 
Cph{T) from the total measured specific heat. The result 
of this calculation is shown in Fig. |^ by the closed circles. 

As mentioned in the introduction it seems quite likely 
that superconductivity in MgB2 is driven by conventional 
electron-phonon coupling. It has been shown that the 
specific heat of conventional, electron-phonon coupling 
driven superconductors is, regardless of the coupling- 
strength, well described by the usual BCS-expression, 



the electronic specific heat at T < Tc may be written 
as 11,0 
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but scaled by the factor as 



A(0) 
1.76fcsTc 



where A(0) denotes the superconducting energy-gap at 
T = K and A is the temperature dependent BCS 
gap-function. This function has been tabulated by 
Miihlschlegel The numerical evaluation of Eq. (H) 

is presented as the dash-dotted line in Fig. |[ We note 
that our data are not well reproduced in this way, but, 
as we show below, the apparent extra contribution is due 
to lattice excitations of the impurities. 

As already mentioned above, our sample contains im- 
purities at concentrations of the order of 1%. Since Boron 
and also Magnesium are rather light atoms, the identi- 
fied impurity atoms are certainly much heavier. Heavy 
impurity atoms may cause resonances in the continuum 
of the phonon excitation spectrum at energies well be- 
low the energy of the lowest optical mode ||T^. In order 
to separate such contributions to the specific heat from 
the electronic specific heat, we fitted our data with an 
expression of the form 

Cp{T)-Cph{T) = Cci{T) + Crcs{T), (6) 

where Cei{T) for T < Tc \s given by Eq. (||) and for 
T > Tc by ^T. To keep the model as simple as pos- 
sible, the contribution to the specific heat due to the 
resonant modes, CresiT), is taken into account by an 
Einstein term similar to Eq. (|^). The impurity concen- 
5|. Therefore, tration ni^p enters as an additional free parameter. The 



fit gives riimp = 0.5%, in line with our expectations, a res- 
onance energy of 30 K, and A(0) = 1.2A:bTc. The ratio 
A(0)/fcsTc is surprisingly small but, as has been shown 
by Swihart jl^ , even for substantially enhanced electron- 

phonon couplings, may be reduced to below the 

original weak-coupling value of 1.76. This is particularly 
the case if low-energy dispersionless phonon modes are 
present. 

The presence of impurities does not significantly alter 
the total specific heat at elevated temperatures. Their 
contribution to the lattice specific heat above Tc is less 
than 10~^ and can safely be neglected in the calculation 
described above in fitting the specific heat at T > Tc- 



ulus i? = 151 GPa and a reasonable value for the 
volume thermal expansion coefficient 6 = 2.5xlO~^K~^ 
at T > Tc, we obtain (Cp - Cv)/T = 1.67 mJ/mole • K^. 
This would lead to a reduction of the electronic specific 
heat coefficient to 7 = 3.83 mJ/mole-K^ and a corre- 
sponding reduction of A to a value of 1.18. According to 
Eq. (^), the ratio A{0) /ksTc would thus be enhanced to 
a value of 1.73. 

The samples used in this work were prepared by S. 
Sigrist. We acknowledge the microprobe analysis inves- 
tigations by P. Wagli. This work was in part finan- 
cially supported by the Schwcizerischc Nationalfonds zur 
Forderung der wissenschaftlichcn Forschung. 



FIG. 4. Electronic specific heat Cei{T) of MgB2 as a func- 
tion of temperature. The solid line represents the rescaled 
BCS-expectation of the electronic specific heat (Eq. (^), and 
the broken line is the electronic specific heat in the hypothet- 
ical normal-state. 



In Fig. Ilwe show the electronic contribution to the spe- 
cific heat in the temperature range between 3 and 75 K, 
which is extracted from the total specific heat by sub- 
tracting the lattice contribution, including the small res- 
onant term discussed above. The solid line in this figure 
is calculated using Eq. (||). We note a rather good agree- 
ment between the calculations and our data, which pro- 
vides further evidence that superconductivity in MgB2 is 
well described by the BCS approximation. 

In conclusion, we have presented a detailed analysis of 
experimental specific heat data for MgB2 , covering an ex- 
tended range of temperatures. The lattice specific heat 
is compatible with a Debye-temperature 9d = 746 K. 
The electronic specific heat is rather well described by 
the BCS approximation, assuming a zero-temperature 
energy-gap of A(0) = 1.2/cbT'c- Finally, we have obtained 
an electron-phonon coupling constant A « 2, which is dis- 
tinctly larger than the usual weak-coupling values. This 
value of A should be considered as an upper limit. If the 
difference Cp — C\/, approximately linear in T, cannot be 
neglected, this would automatically lead to a reduction 
of A. Indeed, with the published value for the bulk mod- 
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